home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DPOSL.z / DPOSL
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPOOOOSSSSLLLL((((3333FFFF))))                                                            DDDDPPPPOOOOSSSSLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPOSL   - DPOSL solves the double precision symmetric positive definite
  10.      system A * X = B using the factors computed by DPOCO or DPOFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE DPOSL(A,LDA,N,B)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAA DOUBLE PRECISION(LDA, N)
  20.         the output from DPOCO or DPOFA.
  21.  
  22.      LLLLDDDDAAAA INTEGER
  23.         the leading dimension of the array  A .
  24.  
  25.      NNNN INTEGER
  26.         the order of the matrix  A .
  27.  
  28.      BBBB DOUBLE PRECISION(N)
  29.         the right hand side vector.  On Return
  30.  
  31.      BBBB the solution vector  X . Error Condition
  32.  
  33.      AAAA division by zero will occur if the input factor contains a zero on the
  34.      diagonal.  Technically this indicates singularity, but it is usually
  35.      caused by improper subroutine arguments.  It will not occur if the
  36.      subroutines are called correctly and  INFO .EQ. 0 .  To compute
  37.      INVERSE(A) * C  where  C  is a matrix with  P  columns
  38.         CALL DPOCO(A,LDA,N,RCOND,Z,INFO)
  39.         IF (RCOND is too small .OR. INFO .NE. 0) GO TO ...
  40.         DO 10 J = 1, P
  41.         CALL DPOSL(A,LDA,N,C(1,J)) 10 CONTINUE LINPACK.  This version dated
  42.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  43.      Subroutines and Functions BLAS DAXPY,DDOT
  44.  
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.